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A comparative study is made on the realizability conditions of 
various special cases as well as the general case of the positive-real 
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I. INTRODUCTION 



The problem of realizability of biquadratic functions as driving 
point immitance has been studied by a number of investigators [1] - [11], 
with realizability conditions given for various special cases as well as 
for the general case of biquadratic functions. It has been proved [10] 
that, given a pair of poles, the location of the zeros are limited by 
two curves, one of which is the locus of the zeros such that the biquad- 
ratic function is a minimum positive real function. 

In Chapter II of this paper, various realizability conditions are 
reviewed and illustrated. FORTRAN programs are written in Chapter III 
for illustration as well as comparison of different cases. 

Many forms of realizations of positive-real biquadratic functions 
are well known. Some of these forms are utilized, in Chapter IV, to 
study the sensitivity problem related to the realizability of the 
function. The method of graphical representation by Chan [13] is applied 
in this study. A computer program is written for the illustration of 
the variation of sensitivity phasors [13] as functions of element values. 

Finally, in Chapter V, a discussion on the results is made together 
with some suggestions for further studies. 
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II. REALIZABILITY CONDITIONS 



A. INTRODUCTION 

Consider a general biquadratic function 



or 



Z(s) 



2 2 
s + 2o s + CO 
z z 



2 2 
s + 2o s + CO 
P P 



(s + z*)(s + z") 
(s + p’) (s + p") 



z(jco) 



2 

-CO 

2 

-CO 



+ 2a ICO + CO 

-zr 

+ 2a ICO + CO 
P 



2 

2 

P 



(co^ - co^) + 

2 2 

(co - CO ) + ico2a 
P P 



( 1 ) 



( 2 ) 



It has been investigated by Chan and Phung [10], that a biquadratic 

function Z(s) as given in equation (1) with conjugate- comp lex zeros in 

the open half plane is realizable as a driving point immitance with 

passive one-port if and only if the zeros are located either on the 

curve C or C. or inside the area delimited by curves C and C., whose 
0 1 0 1 

polar equations are given respectively by 



max CO = CO + 2a cos <p +2 cos (J) (a cos 



z p p 



and 



Min CO = CO + 2a cos , 
z p p z 






z p 



a cos (J) (a cos ^ 
p ^z^ p "z 



* + CO )"] 

z p^ J 

* + 0) )"] 

z p^ J 



(3-a) 



(3-b) 



where (p denotes the compliment of the argument of the zero z , such 
z 

that , 

a = CO cos (p 
z z z 
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B. NECESSARY AND SUFFICIENT CONDITIONS 

To obtain the real part of Z(jcD) of equation (2), multiply and 
divide this equation with the conjugate of the denominator, so that 



Re 




2 2 2 2 2 
(m - o))(a) - oi))+4aa(jj 

_je __z 






^ , 2 
-0) ) + 4o (1) 

P 



(5) 



Z(s) is positive real if the numerator of expression (5) is nonnegative 
for all values of m, that is 



.T/ 2^ 4 2, 2 ^ 2 , .^22 

N(a))=(i) - +0) -4aa)+a)0) 

Z p Z p 2 P 



( 6 ) 



Now we have the following two conditions [10]: 

Condition-1 

2 

N (o) ) has complex or equal positive real roots. 

Condition- 2 
2 

N (o) ) has nonpositive real roots. 

2 

Condition-1 implies the discriminant of N(ai) ) is negative or zero 

i. e. 

D = (o)^ + ( 1 )^ - 4a a “ 4(i)^(u^ ^ 0 (7) 

z p z p z p 
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or after factoring we get 









(o) +o))^-4aa 1(a) - o))^-aa 

zp zpJL^ P ^P 



p ]<0 5 



thus condition -1 is satisfied if 



2 2 
(o) - 0 ) ) 4 : 4a a < (o) +u) ) 

z p z p z p 



(7-a) 



(7-b) 



The right hand inequality is shown as strict since inequality corresponds 
to real roots which would in this case be nonpositive, which is a 
violation of condition- 1 . 

Next, condition-2 implies the following two conditions simultaneously: 
Condition-2a 

D > 0 . 

Condition- 2b 

2 2 

a)+(jj-4aa^ 0 

z p z p 



It can be seen that both conditions 2a and 2b are satisfied if 



4a a X (o) + 0 ) ) 
z p z p 



( 8 ) 



With equations (7-b) and ( 8 ) we obtain the necessary condition for positive 
realness of Z(s) to be 



(o) - 0 ) ) < 4a a 

^ z p^ ^ z p 



(9) 



Now suppose that equation (9) is satisfied with the equality sign. 

2 

Z(s) is then a minimum positive- real function since N(o) ) has a double 
positive root oj^, where 



0). = 0) 0) 

1 z p 



( 10 ) 
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Next suppose equation (9) is satisfied but not necessary as an equality 



then 



a 

z 







4a 

P 



( 11 ) 



The necessary and sufficient condition obtained in the above discuss- 
ion can be conveniently and geometrically interpreted. 

Suppose that the zeros z* and z” are conjungate complex with negative 
real parts and the poles p* and p" are similarily located in the open 
left half plane. 



Let 



a = oj cos 6 
z z 



(12-a) 



where denotes the compliment of the arg\nnent of the zero z*. 



If 0 ) ^ 0 ) , then we let m m + u 

z p z p o 



(12-b) 



If 0 ) ^ oj , then we let m m - u. (12-c) 

z p z p 1 

substituting equations (12-a) and (12-b) into equation (9), and after 
rearranging we get 

2 ' * 

u - 4a (cos d> ) u - 4a a cos (|> <0 (13-a) 

o p ^z o p p z 



and substituting equations (12-a) and (12-c) into equation (9) we get 



u. + 4a (cos ) u. - 4a a cos (|> <0 

1 p z 1 p p z 



(13-b) 



Both polynomials in the left hand sides of the inequalities (13-a) and 
(13-b) have a positive real root and a negative real root, so that these 
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two conditions will be satisfied if and are nonnegative and smaller 
than or equal to the positive-real root of the left hand side polynomials 
in (13-a) and (13-b) respectively. 

Thus solving for the positive real root and then substituting the 

I 

result into equations (12-b) and (12-c) gives, for each value of c() 90°, 

z 

the maximum and minimum values of for which the necessary and suffi- 
cient condition for the positive realness of Z(s) is satisfied. 
Alternatively, we substitute equations (12-b) and (12-c) into equation 
(9) and get 

2 

u - 4o a 0 

z p 

and 

2 

u. - Aa a 0 . 

1 2 p 

Solving these two equations for the positive real root and substituting 

the results into the expression for 03^, gives the maximum and minimum 

values of 03 for which the positive realness of z(s) is satisfied, 
z 

Thus we have in this case the two equations 

max 03 = 03 + 2 /a o (14-a) 

z p z p 

and 

min 03 = 03 - 2 /a a . (14-b) 

z p z p 

Note that these equations will give the same results as equations (13-a) 
and (13-b) when poles of Z(s) are complex, and will also give the 
realizable area of zeros of Z(s) when the poles are real. 

Thus given Z(s) with complex conjugate poles in the left half plane, 
or with negative real poles, Z(s) will be realizable if the zeros of Z(s) 
are complex and lie on or within the curves of equations (14-a) and (14-b) 
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In the special case where Z(s) has real poles and complex zeros, 
it has been shown that z(s) can be realized by a five element series 
parallel circuit [11]. 

The key restriction in this case is given by 

p” (P" - p’) - (z’ - p’')(z" - p") > 0 (15) 

where p', p", z' and z” are poles and zeros of 



Z(s) = (s + z')(s + z") 
' ^ (s + p’) (s + p") 



and p" > p' . 

If we substitute a's and m's for the p's and z's of (15), then we 

get 

I P- ^ 

p = -a + /a - to 

P P P 



fl T 

-a - /a - to 
P P P 



^ p. r 

-a + /a - to 



/2 T 

-a - /a = to 
z z z 



we also obtain 



+ 2o - 0)^ 



2a a 

z p z p 



2 2 ^ . 

u) - to >0 
z p 



rearranging terms, we have 



2 2 

to + 0 ) < 2a a 

z p z p 



+ 2 ^ 
z p p 



Now since 



2 2 2 2 2 
to +(o =(to — to) + 2to to , 
z p z p z p 
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then we get 



2 2 
(o) - 0 ) ) ^ 2a a + 2a /a - 0 ) - 2o) w 

zp zp zp p zp 



Comparing this condition to that for the general case in equation (9) , 
we obtain 



so that this is indeed more restrictive than the general case and we 
should not be surprised to find the area of realizability for this case 
enclosed by that of the general case. 

In fact this argument could be used as a proof of the realizability 
of Z(s) given in (15) , since the area of this case is always enclosed 
by that of the general case as shown above. 

The construction of the realizability region in this case will be in 
accordance with reference [11] , that is 



2a a 

z p z P P z p 





2 

(a> - a> ) < 2a a + 2a 

z p' z p z p p 




> 0 ) < 4a a 

z p z p 



i i 

pM-r-p (p-p) , 



where p" > p’. 



C. SAMPLE ILLUSTRATIONS 



(a) Complex poles > a^) 



Let a = 3 
P 



and 0 ) = 5 . 

P 



14 



Then from equations (14-a) and (14-b) we get 



and 



(jj max = m + 2/a a 
z p z p 



= 5 + 2/3^" = 5 + /12a , 

z z 



(D min - 0) - 2/a a 

z p z p 



5 - /12a . 

z 



We obtain the limitation of a as follows: 

z 



a^ < for complex zeros. 



If we let 0) max = a , then we get 
z z 



u = 5 + /12a 
z z 



= 20.8 = a , maximum value for a , for t 
z’ z 



Similarly if min = a^, then 



0) = 1.2 = a , maximum value for a , for m 

z z z : 



(b) Real Poles 



Let 



a =2.5 and 
P 



= 2 with poles at (-1, -4). 



Then we obtain 



(1) max = 0) + 2/ a a 

z p z p 



I max. 
z 



min. 



= 2 + /10a 







I 



uj. 







and 



iii min * 0 ) - vlOa 

z p z 

If 0 ) max = a , then 
z z* 

(1) = 2 + /lOa 

z z 

= 13.7 = a^, max for o)^ max and complex zeros. 

From the examples given in this section it can be seen that the 
realizability region of a biquadratic function is dependent upon the 
restriction we place on the function. That is, as we place more 
restrictions on the function and proceed from a general function with 
complex poles to one with real poles which is realizable by certain 
technique, we see that at each step the realizability region decreases. 
This decrease of freedom from the additional restriction is to be 
expected. 
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III. COMPUTER STUDY ON REALIZABILITY 



For a given location of the complex pole pair (a and w specified) 

f 

in the open left half plane, and are simply plots of versus 

Thus the computer program was written to accept and as inputs and 

to generate and on either the line printer or the Calcomp plotter. 

Some sample outputs are shown in Figure A-l.A and A-l.B in Appendix A. 

The curves are for a normilized set of pole locations with 

0 ) = 1 
P 

and 

0° < ((.' < 90°, 

P 

? T 

where <b is the compliment of the argument of the pole p . 

P 

Varying (f)^ from 0^ to 90° corresponds to shifting the poles along 
a circular path from the negative real axis to the jm axis. 

The program was then modified to produce a set of curves and 

for various pole locations. Given a specific pair of poles, two types 
of plots may be obtained, depending upon the requirements of the problem. 

If an accurate plot of the realizability region is required, the 
program can be used so that a plot of the realizability region is 
required, the program can be used so that a plot similar to that shown 
in Figure A-l.A is obtained. Otherwise if a quick estimate is all 
that is desired, a plot similar to the one in Figure A-l.B can be used. 
Other outputs were obtained for various pole locations. 

Results obtained for the following cases are included, with 
program listing, in Appendix A and B respectively; 

1. Complex poles with moves toward a constant 



17 








i 



I 



■ii 



2. Complex poles with moves toward a constant 

3. Complex poles that are moving toward the origin. 

4. Real negative poles with one of which moving toward the origin 
while the other remains stationary. 

5. Real negative poles with one of which moving toward infinity 
while the other remains stationary. 

6. Double real poles moving from the neighborhood of the origin 
toward infinity. 

7. Pure imaginary poles moving toward the origin. 

8. Pure imaginary poles moving away from the origin. 

The above cases are now described in the following sections. In 
each case a figure is included showing the pole locations with their 
values shown in the accompanying table, in which the location of the 
corresponding computer output plots are indicated. 
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Case 1 



This case illustrates the zeros region for given complex poles with 
a moves toward a constant w . 




Five different poles ( conjugate poles are not shown ) are given in 
accordance to show the zero regions for each as tjp remains constant. 
The locations of each pole are tabulated in Table 1. 



Table 1 



Pole 




(13 

P 


Computer Plot 




0.999 


1.0 


Figure A-1.A,B 


^2 


0.707 


1.0 


Figure A-2 




0.55 


1.0 


Figure A-3 




0.25 


1.0 


Figure A-4 




0.1 


1.0 


Figure A-5 
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Case 2 




This case illustrates the zeros region for given complex poles with 
u)p moves toward a constant • 

Five different pole locations are shown in Table 2, 



Table 2 



Pole 






Computer Plot 


C 

1 


1.0 


1.05 


All plots are 


^2 


1.0 


2.0 


shown in 


C 

3 


1.0 


3.0 


figure A-7 . 


C 

4 


1.0 


5.0 






1.0 


8.0 





20 



I 



I 



Case 3 



\ 





' ( ju>) 


;\ 

S ; 




S ; : 




C 4 I ; , ♦ 

' 1 1 \ 




"^5 ; : i : \ 


CJ3 


^Pa', 

''P4' 





^P4\ 

Figure 4. 



This case illustrates the zeros region for given complex poles that 
are moving toward the origin. Five different pole locations are 
shown in Table 3. 



Table 3 



Pole 


T 

P 


0 ) 

P 


Computer Plot 




4 .0 


30.652 


All plots are 




2.0 


15.326 


shown in 




0.5 


3.834 


Figure A- 8 . 




0.1 


0.766 






0.01 


0.077 
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Case 4 



This case illustrates the zeros region for given real negative 
with one of which moving toward the origin while the other remains 
t ionary. 




Figure 5. 



Five different pole locations are shown in Table 4. 



Table 4 



Condition 


Pole 




Pole P^ 
2 








8.0 


, 0 


4.0, 0 


6.0 


5.657 


s 


8 .0 


, 0 


1.0, 0 


4.5 


2.828 


S 


o 

00 


, 0 


0.5, 0 


4 .25 


2.0 




8 .0 


, 0 


0.1, 0 


4.05 


1.414 




8 .0 


, 0 


0.05,0 


4.025 


0.632 



All plots are shown in computer output Figure A-9. 



poles 

sta- 



■5S 
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Case 5 



This case illustrates the zeros region for given real negative poles 
with one of which moving toward infinity while the other remains stati- 
onary. 



/fjtj 







P 



2 



<r 

p 



Figure 6 . 



Five different pole locations are shown in Table 5 . 



Table 5 



Condition 


Pole 


Pole 


P 

2 


(T 

P 


‘^P 


C 


0 

0 

00 

■ 


2.0, 


0 


5.0 


4.0 


1 


C 


15.0,0 


2.0, 


0 


8.5 


5.477 


2 


C 


30.0,0 


2.0, 


0 


16.0 


• 7.746 


3 


"4 


50.0,0 


2.0, 


0 


26.0 


10.0 


C 


100.0,0 


2.0, 


0 


56.0 


14.142 



All plots are shown in computer output Figure A-10. 
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Figure 7. 



This case illustrates the zeros region for given double real poles 
moving from the neighbourhood of the origin toward infinity. Five 
different pole locations are shown in Table 6. 



Table 6 



Condition 


Pole P^-P 2 


1 


u) 

P 


Computer Plot 




1.0, 0 


1.0 


1.0 


All plots are 


^2 


2.0, 0 


2.0 


22.0 


shown in compu- 




5.0, 0 


5.0 


5.0 


ter output 




10.0, 0 


10.0 


10.0 






25.0, 0 


25.0 


25.0 
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Case 7 




Figure 8. 

This case illustrates the zeros region for given pure imaginary poles 
moving toward the origin. Five different pole locations are shown in 
Table 7. 

Table 7 



Condition 


<r 

p 


cO 

P 


Computer Plot 


C 

1 


0.0 


3.0 


All plots are 


<^2 


0.0 


1.0 


shown in 


S 


0.0 


0.5 


Figure A-12 . 


"4 


o 

o 


0.1 






0.0 


0.05 
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Case 8 




This case illustrates the zeros region for given pure imaginary poles 
moving away from the origin. Five different pole locations are shown in 
Table 8. 



Table 8 



Condition 


a 

p 




Computer Plot 




0.0 


4.0 


All plots are 


<=2 


0.0 


10.0 


shown in com- 




o 

o 


<.C . 0 


puter output 




0.0 


100.0 


Figure A- 13. 




0.0 


200.0 
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IV. SENSITIVITY STUDY 



A. INTRODUCTION 

It has been shown [12], that in a linear and not necessarily 
reciprocal network, with only one independent source and initially at 
rest, the ratio of the Laplace transform of the current in any branch 
or of the voltage across any node pair to the Laplace transform of 
either an independent voltage or current source is given by 
W T (0) + X T («) 

"W = - - w^V - (1^) 



where 



X 

T(x) 



T (0) 

X 



T (») 

X 



an adjustable immitance. 

the network function relating a response to an 
exitation. 

Lim T(x) = the network function evaluated with 

X — *■ o 

X = 0, is assumed to be finite. 

Lim T(x) = the network function evaluated when 

X — ► ® 



X = ®, is assumed to be finite. 

W = the Thevenin immitance (independent of x) seen looking 
back into network from the terminals of the adjustable 
immitance. W is assumed to be finite. 

Equation (16) is generic in form in that x can be either an impedance 
or admittance. In this paper the symbol T(x) is used to denote a driving 
point immitance. 
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B. GRAPHICAL REPRESENTATION OF T(x) 

Consider equation (16) for the case of a driving point immitance 
function as shown in Figure 10. 

Dividing numerator and denominator of equation (16) by W, we get 



T(x) 



X 

T (0) + T (oo) 

w ^ 



X 

1 + 

W 



which can be rewritten as. 



T(x) 



X 

T^(oo) (1 + ) + T^(0) - T^(co) 



X 

1 + — 
W 



or 

T (0) - T (->) 

T(x) = T^(co) ^ 

X 

1 + 

W 



(17) 



Let 



T(x) = T^(«) + V 



(18) 



where 



V 



T (0) - T («) 

X X 

X 

1 + — 

W 



(19) 
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\ 




T(x). 




Figure 10. 

Driving p6int immitance with variable 
component x. 
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From equation (19) we have 

X 

V + V = T (0) - T (co) (20) 

y X X 

we note that ^nd ^.re fixed quantities. 

Now treat each quantity in (20) as phasor, i.e. 

V = |v| exp. 

and 

W - |w| exp. j 62 

then the left hand side of equation ( 20 ) represents the sum of two phasors. 

If X is real, then the phasor 
X V 

V = X I 1 exp. j(0 - 9 ) 

W W ^ ^ 

is lagging behind the phasor V by an angle 
9 = 8 ^- 9 ^. 

Since fixed quantities, the right hand side 

of equation ( 20 ) is constant so that the left hand side of ( 20 ) is 
constant. This requires that the angle must be constant (as x varies). 

This implies that the tip of phasor V describes the arc of the 
circle, of which the quantity ^ chord as shown in 

Figure 11. 

To construct the circle we proceed as follows [13] : 

1. Plot phasors from the origin. 

2. Draw a straight line through the tips of the two phasors in 
step 1 . 
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3. Calculate the angle and draw a straight line from the tip of 
T^(0) forming angle with having the value -0 measuring clock- 



Therefore is tangent to the circular locus of T^. 

The centre of this circle, 0*, is located at the intersection of 
the two perpendicular lines erected on the chord in the middle of 
T^(0) - at the tip of shown in Figure 12. 

Knowing the locus of T(x) as x varies from zero to infinity, we 
can graphically determine T(x) at any value of x, 

C. GRAPHICAL REPRESENTATION OF SENSITIVITY PHASOR 
If we define the sensitivity as 

T(x) ^ d T(x) 

X dx 

then with equation (16) we can obtain, [13] 



wise from D. 



1 



gT(x) 



T^(0) - 



^2 

X 



W 1 + — 



W 



or 




V 



X 

V 

W 



T^(0) - T^(«) 



where 



V 



Tx(0) - T^(«) 



X 



1 + 

W 
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Figure 12- 



Graphical representation of locus of T(x), 




Graphical representation of sensitivity 
phasor S, 
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It can be shown [13] that the sensitivity of a function may be 
graphically represented by a phasor such as PS as shown in Figure 13. 

D. SAMPLE ILLUSTRATION 

As an illustration of the above discussion in this chapter, let us 
consider a biquadratic p.r. function as an RC driving point admittance, 
which can be realized in Foster 2 form as shown in Figure 14. 

Let 



Y(s) 



2 

s + 4s + 3 
s^ + 6 s + 8 



(s + 1) (s + 3) 
(s + 2) (s + 4) 



which can be realized in Foster 2 form with the element value as given 
in Figure 5, in which all resistance and capacitances are measured in 
Ohms and Farads respectively. 

Now let R 2 change by the amount AR, such that 

0 < AR < w 

Thus, Figure 15 is modified as shown in Fgiure 16. 

We can calculate T^(0) , W as follows: 



T^(0) = Y(s, AR=0) 




+ 



1 

+ 1/sC^ 



+ 



1 

R^ + l/sC^ 



07 c 4 . «125 s . .09375 s 

.5 s + 1 .25 s + 1 



s^ + 4s + 3 
s^ + 6 s + 8 



. (24) 



T^(«) = Y(s, AR=») = — + r - Vi/ T c -- 

oil 



= .375 + 



.125 s 
.5 s + 1 



2.5 s + 3 
4 s + 8 



(25) 
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Figure 14. 

The realization of Foster-2 form of 
driving point admittance Y(s). 




Figure 15. 

The realization of biquadratic p.r. driving point 
admittance in Foster 2 form, where 

Y(s) = _ 

s^+6s+8 
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sC^ Y(s, AR=co) 



. . 10.667 . 4 s + 8 

2.667 + ’s + 3 



14.5 s + 40 s + 99 

7.5 + 9 s 



( 26 ) 



From expression (26) we obtain the expression to be 



X 

W 



"1.5 s X + 9 s X 



(27) 



14.5 s + 40 s + 99 
Substitute equations (24), (25), (26) and (27) into equation (22) we get. 



T(x) 



s + 4s + 3 2.5 s + 3 

S^ + 6s + 8 4 s + 8 



14 



.5 + 40 s + 99 ^ 7.5 xs^ + 9 xs ~ |^ 

7.5 + 9 s L 14.5 + 40 s + 9?J 



, (28) 



then for simplification we write equation (28) above as. 



T(x) 

S 

X 



N(s,x) 

D(s,x) 



(29) 



where N9s,x) and D(s,x) represent the numerator and denominator of the 
equation (28) respectively. 

In the frequency domain, the expression in equation (29) can be 
written as , 



T(x) 

S (jw,x) = 

X 



N(.1cl),x) 

D(jw,x) 



(30) 
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Figure 16. 

Driving point admittance Y(s) with variable resis- 
tance 4R, 



1 w 


LdRi 

Li.. 


1 

1 

1 

1 

Y(s) ; ^ 


'XR 

o 




1 i— 1 

J i 1 


1 

1 

1 

1 

6 


1 











Figure 17. 

Driving point admittance Y(s) with variable resis- 
tance R^' = ^2 ^ 
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= _ Re lN(jg>,x)] + j Im [N(.ito,x)] 

Re [D(j(i>,x)] + j Im [D(ja>,x)] 

where Re [N(jw,x)], Re [D(jo),x)], Im [N(jw,x)] and Im [D(jm,x)] are 
real and imaginary parts of equation (30) repectively. 

From equation (30a) we can obtain the magnitude of that sensitivity 
function. If we write equation (30a) in exponential form we have 



sT(x) 

X 



(jw,x) 



is!“i 



exp 






where 



sKx) 

X 



magnitude of equation (30a) 



{Re [N(jo),x) ] + {lm[N(j(i),x) ] }^ 

(Re [D(jco,x)] ^ + {Im[D(j(i),x) ] }^ 






(31) 



= the argument of the numerator of equation (30a) . 

= the argument of the demominator of equation (30a) . 

E. COMPUTER STUDY ON SENSITIVITY 

In this section, the sensitivity of the function T(x) with respect to 
the variation of one of its components will be investigated by computer. 
Consider the general expression Y(x) in equation (16), which is repeated 
here for convenience. 

W T (s,x = 0) + X T (s, X = «>) 

T(x) = . (16) 

W + s 



The sensitivity function of this equation is 



eT(x) _ 



T^(s, X = 0) 



T (s, X = 
x' ’ 



« f + -f-] 



( 22 ) 
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As an example let us consider the function Y(s), as we derived in section C 



Y(s) 



(s + DCs + 3) 
(s + 2) (s + 4) 



and the realization of this function in Foster-2 form is shown in Figure 
17. 

For our investigation, let change with the amount of AR in series 
with R^ and let 

R^' * ^2 

T (s R * ) 

To obtain the sensitivity S * 2 in terms of frequency and R«’, 

^2 ^ 
we use the same calculation from equation (30) . 

The computer output in Appendix B shows us the plot of the magnitude 

of sensitivity function (30) versus ^ 2 '* this purpose we use equation 

(31) as reference. 

Various cases are investigated using computer, resulting in plots 
of the magnitude of equation (31) versus R^’ for different frequencies. 
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V. CONCLUSION 



In this paper, a comparative study has been made on the realizability 
conditions of a general biquadratic function. The graphical method of 
Chan and Phung [10] was used for this study and computer programs are 
written for verifications. General observations can be made as follows: 

1. The realizability for the zeros expands considerably as the poles 
move from the immaginary axis to the real axis. This results both from 
C^ moving away from the origin and C^ moving toward the origin. 

2. The limits on C^ and C^ are found easily from equations (3a) and 

(3b) . As Op approaches zero, corresponding to the poles moving toward 

the 10 ) axis, both C and C. tend toward the circle o) = o) for all (p \ 

As can be seen from Figure A-6, C^ and C^ separate very rapidly for 

small movements of the poles away from the jm axis. We can conclude that 

when a = 0 we require o) = o) in order for Z(s) to be realizable with 
P ^ z p 

complex zeros. This can also be seen easily from equation (9). 

3. As Op approaches u)^, corresponding to the poles moving toward 

the real axis, C approaches a value of -5.828 (d on the real axis, and 
o p 

Cj^ similarly approaches a value of -0.172o)p, as we can see in Figure A-l.A. 

A brief study was made on sensitivity functions [13], some computer 
results are obtained. However, due to time limitation, a complete 
investigation in sensitivity was not possible and is left as a 

suggested topic for further studies. 
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APPENDIX A 



COMPUTER OUTPUT 



A. ZERO’S REGION 

The illustration of zeros regions will be shown in this appendix. 

The computer outputs show us the regions for various given omega 
and sigma poles. 

Various cases are explained as follows: 

Figure A- la and A- lb show the zeros regions for given sigma pole 
= 0.999 and omega pole = 1.0 using Calcomp and PLOTP in the computer 
program respectively. 

All the following plots show us the computer outputs by using Calcom 
in the computer program. 

Figure A-2, A-3, A-4 and A-5 show us the plots of zeros regions for 
given sigma poles equal to 0.707, 0.55, 0.25 and 0.1 respectively with 
the amount of omega pole equal to 1.0 each. 

Figure A-6 shows all plots given in Figure A-la, A-2, A-3, A-4 and 
A-5 in one computer output. 

Figure A-7 shows the plots of zeros regions as pole moves at sigma 
pole constant. 

Figure A- 8 shows the zeros regions as pole moves to the origin. 

Figure A-9 and A-10 show the plots of zeros regions as pole moves 
to the origin and to infinity respectively. 

Figure A-11 shows the plot of zeros regions of a double real negative 
poles as both move along the negative axis. 

Figure A-12 and A- 13 show us the zeros region plots of a given pair 
of pure negative poles, as both move to the origin and to the infinity 
respectively. 
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B. SENSITIVITY FUNCTION 



Various cases of sensitivity as a function of are illustrated in 
the computer outputs for different particular frequencies. 

The computer outputs show us the plots of sensitivity versus R* ^ 
and the tabulations of this plot. 

The resistance R ’2 is select to be varied from the amount of 1 Ohm 
to 50 Ohms with increment to 1 Ohm. And the frequencies are choosen 
beginning at 1 rad/sec. up to 10,000 rad/sec. 

Figure B-1, B-2, B-9 show us the sensitivity functions as 

R *2 moves from 1 to 50 Ohms, at the frequencies of 1, 5, 10, 50, 100, 
500, 1000, 5000 and 10000 rad/sec. The tabulations of those plots 
are also shown in Table B-1, B-2 B-9 respectively. 

Figure B-10 shows us five different plots of sensitivity functions 
illustrated in one computer output. 
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A-1 B 



This plot shows the Zero Eegion for given Sigma pole 
« 0.999 and Omega pole » 1.000. 
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The plot of eensltivity va R^' la shown for given 
frequency 1 r-f’/aec. 
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Sensitivity plot vs R^' is shown here for given 
frequency 5 rad/sec. 
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for given frequency 10 rad/sec. 
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The plot of eensltlvlty ve R^' is ehown here for gi' 
frequency 50 rad/eec. 
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plot of eenBitlvity vs R^' for given frequency 
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frequency ^0 rad/eec. 
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frequency 1000 rad/sec, 
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The plot of oeneltlvlty vs R. 
given frequency 5000 rad/sec, 
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The plot of Bensltlvit^ vs R2' 
given freouency 10000 rad/sec. 
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plots of sensitivity vs R. 
given various frequencies, 



The V4LUE CF SENilTIVITT i&T GIVEN W = 



P2S 


SENSITIVITY 


C . 0 


-0.005995 


l.GCO 


-0.005160 


2.000 


-C.00AA89 


2. COO 


-0.0C39A0 


A. COO 


-0.003A36 


5.0C0 


-0.002107 


6. COO 


-0.002736 


7.000 


-0.002512 


8.000 


-0.002277 


c.OOO 


-C. 002073 


1C. COO 


-0.001896 


11.000 


-0.0017A0 


12. COO 


-C.C01603 


12.000 


-0.C01A81 


lA.CCC 


-0.001373 


15.CCC 


-0.001276 


16. COO 


-0.001189 


17.000 


-0.001111 


18. COO 


-C.OCIOAO 


1^.000 


-O.OOC976 


2C.CC0 


-0.000917 


21.000 


-C.0C086A 


22.000 


-0.000815 


22.CC0 


-C.00C770 


2^. COO 


-C. 000729 


25.000 


-C.00C691 


26.CC0 


-C.000656 


27.000 


-C.000622 


28.000 


-0.000593 


2S.CCC 


-C.C00565 


2C.C00 


-0.000539 


21.000 


-0.000515 


22.CC0 


-C.GC0A92 


22.000 


-C.OOOA71 


2^.CC0 


-0.000A51 


25. COO 


-C.000A32 


26.000 


-0.000A15 


27.000 


-0.000398 


38.000 


-C.000383 


3S.CC0 


-0.000268 


AC. COO 


-C.00025A 


Al.OOC 


-0.C002A1 


A2.000 


-0 .000229 


A3. COO 


-C. 000217 


AA.COO 


-0.000206 


A5.000 


-0.000296 


A6.000 


-C. 000286 


A7.0CC 


-0.000276 


A8.CC0 


-0.000267 


AG. COO 


-0.000259 


5C.0C0 


-0.000251 



TABLE B 



1.000 
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THE VALUE CF SENEITIVITV AT GIVEN W 



5,00C 



R2S 


SENSITIVITY 


G.O 


-0.086060 


l.CCO 


-0.054744 


2.CC0 


-0.037862 


3,CCO 


-0.027736 


4.000 


-0.021188 


5.0CC 


-0.016 711 


6. COO 


-0.013516 


7.000 


-0.011156 


e.oco 


-0.009365 


G.COO 


-0.007972 


IC.QCO 


-0.006868 


ll.CCC 


-0.005979 


12.000 


-0.005252 


12.000 


-0.004649 


14.000 


-0.004145 


15.000 


-0.003718 


16.0CC 


-0.003354 


17. COO 


-0.003041 


18. COO 


-0.002770 


IS.CCC 


-C .002534 


2C.0CC 


-C. 002326 


21.000 


-0.002143 


22.000 


-0.C01981 


23.000 


-C .001837 


24.000 


-C. 001707 


25.000 


-0.001591 


26.000 


-C. 001487 


27.000 


-0.001392 


28.000 


-0.001306 


29.000 


-C. 001228 


30.000 


-0.001157 


31.000 


-C .001091 


32.000 


-0.001031 


33.000 


-C. 000976 


34.000 


-C. 000925 


35.000 


-C.000878 


36.000 


-C.000835 


37.000 


-C. 000795 


38.000 


-0.000757 


39.000 


-0.000722 


40.000 


-C.00C690 


41.000 


-0.000659 


42.0CC 


-C. 000631 


43.000 


-C.C00604 


44.000 


-C. 000579 


45.000 


-0 .000556 


46.000 


-C.C00534 


47.000 


-C. 000513 


48.000 


-C.C0C494 


49.000 


-C.G00475 


5C.CCC 


-C. 000458 



TABLE B ->2. 



67 



THE V^LUE CF SENSniVITT /ST GIVEN W 



10.000 



P2S 


SENSITI VITY 


0 • c 


-C. 227638 


l.CQO 


-'C.122C24 


2.0C0 


-0.076772 


2.CC0 


-0.052258 


^•CQO 


-0.037941 


5.CQQ 


-0.028735 


6. COO 


-0.022506 


7. COO 


-0.018097 


8.0CQ 


-0.014865 


«.CQC 


-0.012425 


IC.OCO 


-0.010540 


ll.COQ 


-0.009052 


12.0CC 


-C.0C7858 


13.000 


-0.006885 


lA.CCC 


-C. 006082 


15.000 


-0.005412 


16. COO 


-C .004846 


17.000 


-C. 004365 


18.0C0 


-C. 002952 


IS. 000 


-0.003595 


20.CCC 


-0.003284 


21.000 


-0.003012 


22.000 


-0.002772 


22.000 


-0.002560 


2^.000 


-C. 002271 


25.000 


-0.C02202 


26.CC0 


-C.002C51 


27.000 


-0.001915 


28.000 


-0.001792 


2S.000 


-0.0C1680 


20. COO 


-0.001579 


21.000 


-0.001486 


22.CC0 


-O.C01401 


23.000 


-C .001224 


24.000 


-0.001252 


25.0C0 


-0.001187 


26. COO 


-0.001126 


27.0C0 


-0.001070 


28.CC0 


-0.001018 


2S.CC0 


-C. 000969 


4C.0Q0 


-0.000924 


41.000 


-0.000882 


42.0C0 


-0.000843 


42. COO 


-C.C00807 


44.000 


-0.000772 


45.CC0 


-C.C0C740 


46.CC0 


-C.CCG71Q 


47.000 


-C .000682 


48. COO 


-0.000655 


4S.CC0 


-0.000620 


5C.CCC 


-0.000606 



TABLE B - 3 
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THH VALUE CF SENSITIVITY AT GIVEN W = 50.000 



P2S 


SENSITIVITY 


c • c 


-C.67A758 


1.000 


-0.270A19 


2.000 


-0.1AA071 


3.000 


-0.C89189 


4.CC0 


-C. 060562 


5.000 


-C.0A3783 


6. COO 


-0.032117 


7.000 


-C.025S19 


8.CC0 


-0.020626 


S.COO 


-0.017112 


IC.OCO 


-0.01A3C5 


ll.CCO 


-0.012135 


12.000 


-C.010A2A 


13.000 


-0.009050 


14.CC0 


-0.007932 


15.0CC 


-0.CC7CC8 


16. COO 


-0.006237 


17.000 


-0.005587 


le.coo 


-C.C05023 


IS. 000 


-0.00A557 


2C.CC0 


-0.00A1A6 


21.000 


-0.003788 


22.CC0 


-0.003A75 


23.000 


-0.003199 


24.CC0 


-0 .C0295A 


25. COG 


-C. 002737 


26. COO 


-0.0025A2 


27.000 


-0.002368 


ze.ccc 


-C. 002211 


2S.0C0 


-C. 002069 


3C.0C0 


-0.0019A1 


31.CC0 


-C.C0162A 


32. COO 


-C. 001717 


33.000 


-0.001619 


3A.CC0 


-0.001530 


35.000 


-0.001AA7 


36.CC0 


-0.001371 


37.CC0 


-0.001301 


36.CCC 


-C. 001227 


3S.CCC 


-0.001176 


AO. COO 


-0.001121 


Al.OCC 


-0.001069 


A 2 . 0 CO 


-0.001020 


A3.CC0 


-C. 000975 


AA.OOO 


-0.000922 


A5.C00 


-C.C00893 


A6.000 


-0.000856 


A7.0C0 


-0.000821 


A8.000 


-0.C00789 


AS.CCC 


-0.000758 


5C.CCC 


-0.000729 



TABLE B - 4 
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THE VALUE CF SENSITIVIIV AT GIVE^ W 



100,000 



R2S 


SENSITIVITY 


0,0 


-0.739854 


1.000 


-0.284703 


2.0CC 


-0.149252 


3.CC0 


-0.091733 


A. COO 


-0.061999 


5.00C 


-0.044683 


6.CC0 


-0,022724 


7.CC0 


-0.0 26252 


8.CC0 


-0.021157 


c.CCO 


-0,017259 


10.000 


-0.014499 


11.000 


-0.012292 


12.CC0 


-C.01C553 


12,000 


-0.009158 


l^.CCO 


-C.C08C23 


15.000 


-0.C07C66 


16.000 


-0,006204 


17.000 


-0,005645 


lE.QCO 


-0.005084 


l^.OCC 


-0.G04603 


2C.C00 


-0.004187 


21.000 


-0.002625 


22.000 


-0.003508 


23.000 


-C. 002229 


24 .000 


-0,002981 


25.000 


-0,002762 


26.000 


-0.002565 


27.000 


-0.002389 


26.CCC 


-0,002230 


2S.C00 


-C.C02C87 


20.000 


-0.001957 


21.000 


-C.G01839 


22.000 


-0.001731 


22.000 


-0.001633 


24. COO 


-0,001542 


35.000 


-0.001459 


26.000 


-C.C01283 


37.000 


-C. 001212 


28.000 


-0.001247 


2^.000 


-C. 001186 


40.000 


-C.COl 120 


41.000 


-0.001077 


42,000 


-0.001028 


42. COO 


-C.0C0983 


44,000 


-O.OC0940 


45.000 


-C.CC0900 


46.000 


-0 .000862 


47.000 


-C. 000828 


48.CCC 


-C.C0C795 


4C.0C0 


-0.00C764 


50. COO 


-C. 000734 



TABLE B ~ 3 
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THE VALLE CF SENSITIVITV AT GIVEN M = 500.000 



R2S 


SENSITIVITY 


0.0 


-0.764546 


l.CCO 


-0.289726 


2.0CO 


-0.151158 


2. GOO 


-0.092591 


4. COO 


-0. 062480 


5.0C0 


-0.044982 


6. COO 


-0.032925 


7. COO 


-0.026495 


8.000 


-0.021263 


^.CCC 


-0.017441 


IC.OCC 


-0.014563 


ll.CCO 


-0.012244 


12.000 


-0.010595 


13.CCC 


-0.C09193 


14.0C0 


-0.008053 


15.CCC 


-0.007112 


16.C0C 


-0.006227 


17.0C0 


-0.005665 


16.000 


-0.005101 


19.CC0 


-0 .004618 


2C.CC0 


-0.004200 


21.0CC 


-0.003637 


22. COO 


-0.003519 


23.000 


-0.002238 


24.000 


-C.0C2990 


25.CC0 


-O.C02770 


26.000 


-0.002573 


27.000 


-0.002296 


26.000 


-0. 002237 


2^.000 


-0.002093 


50. COO 


-0.001963 


31.000 


-0.001644 


32.000 


-0.001736 


22.000 


-0.001637 


24.000 


-0.001546 


35.000 


-0.001463 


26.000 


-0.001286 


27.000 


-0.001215 


38.000 


-0.001250 


3S.C00 


-0.001189 


40. COO 


-0.001132 


41.000 


-0.001080 


42.000 


-0.001031 


42.000 


-0.000985 


44.000 


-0.000942 


45.000 


-0.000902 


46.000 


-C.OOG865 


47.000 


-C .000620 


46.000 


-0.000797 


49.000 


-0.000765 


50.000 


-0.000726 
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THE VALUE CF SENSITIVITY AT GIVEN W = 1000.000 



P2S 


SENSITIVITY 


€.0 


-0.765354 


l.OCO 


-0.289887 


2.CC0 


-0.151215 


3.CC0 


-0.092619 


^•COO 


-0.062495 


5.C0C 


-0.044992 


6.000 


-0.033921 


7.QC0 


-0.026499 


6.CCC 


-0.021266 


S.OCO 


-0.017^43 


1C. 000 


-0.014565 


11. COO 


-C.C12345 


12.000 


-0.010596 


13.CQC 


-0.009195 


l^.OCO 


-C.008C54 


15.0CC 


-C.C07113 


16.0C0 


-0.006327 


17. COO 


-0.005665 


18.000 


-0.005102 


IS.CCO 


-0.004619 


2C.CCC 


-C.CG4201 


21.000 


-0.003827 


22.CCC 


-0.003519 


23.CC0 


-0.003239 


2^.CCC 


-0.002991 


25.CCC 


-C.CC2770 


26.CCC 


-0.002573 


27.000 


-0.002396 


28.000 


-0.002237 


2G.C00 


-0.002093 


3C.CCC 


-0.001963 


31.000 


-0.001844 


32.000 


-0.001736 


33.CCC 


-0.001637 


34.C0C 


-0.001546 


55.CC0 


-0.001463 


36. COO 


-0.001386 


37.000 


-0.001315 


38.000 


-0.001250 


39.000 


-0.001189 


40.000 


-0.001133 


41.000 


-0.001080 


42.000 


-0.001031 


43.000 


-C. 000985 


44.000 


-0.000943 


45. COO 


-0.000903 


46.000 


-0.CO0865 


47.000 


-C. 000830 


48.000 


-0.0CC797 


49.000 


-0.000765 


5C.CCC 


-0. COO 736 
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THE VALUE CF SENSITIVUV AT GIV£^ W = 5000.000 



P2S 


SENSITIVITY 


c.c 


-0.765613 


l.COC 


-0.289929 


2.0GC 


-0.151224 


2.00C 


-0.092627 


4.CC0 


-C. 062500 


5.CCC 


-0.044995 


6.CCC 


-0.033923 


7.CC0 


-C. 026501 


8.000 


-0.021267 


9.QC0 


-0.017444 


IC.CCC 


-0.014566 


11.000 


-0.012246 


12.0CC 


-0.010597 


12.000 


-C. 009 195 


14.0CC 


-0.008054 


15.000 


-0.007113 


16.CCC 


-0.C06327 


17.0CC 


-0.005665 


16.000 


-0.005102 


Ic.CCO 


-0.004619 


2C.CCC 


-0.004201 


21.CC0 


-0. 003627 


22.000 


-0.002519 


23.000 


-0.002239 


24.000 


-0.002991 


25.000 


-0.002770 


26.CCC 


-0.002573 


27.000 


-0.002396 


28.000 


-0.002237 


29.000 


-C. 002093 


20. COO 


-0.001963 


21.CC0 


-0.001844 


22.000 


-0.001736 


33.000 


-0.001637 


24.000 


-0.001547 


25.000 


-0.001463 


26.000 


-0.001286 


37.000 


-0.001315 


26.CC0 


-0.001250 


39.000 


-0.001189 


40.000 


-0.001133 


41.000 


-c.ocioao 


42.000 


-0.001021 


42.000 


-0.000985 


44.000 


-0. 000942 


45.000 


-C. 000903 


46.000 


-0.000665 


47.000 


-C.CC083C 


46.000 


-0.C0C797 


49.000 


-0.000766 


50.000 


-0.000736 
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THE VALUE CF SENSITIVm AT GIVEN W = 10000.000 



R2S 


SENSITIVITY 


C • 0 


-0.765621 


I.OOO 


-0.289940 


2.0CC 


-0.151234 


3.CC0 


-0.092627 


4.0C0 


-0.062500 


5.CC0 


-0.044995 


6.CG0 


-0.033933 


T.OCO 


-0.026501 


8.000 


-0.C21267 


<:.CCO 


-0.017444 


IC.CCC 


-0.014566 


11.000 


-0.012346 


12.000 


-C. 010597 


13.CC0 


-0.009195 


14.000 


-0.008054 


15. COG 


-O.OC7113 


16.CCC 


-0.006327 


17.0C0 


-0.C05665 


18.CC0 


-O.C05102 


l^.CCO 


-0.CC4619 


2C.CCC 


-C.C04201 


21.000 


-C .003837 


22.000 


-0.003519 


23. COO 


-0.C03239 


24.CCC 


-0.002991 


25.000 


-C. 002770 


26. COO 


-0.002573 


27. COG 


-0.002396 


28. COG 


-C.C02237 


29. COO 


-0.002093 


3C.CCC 


-0.001963 


31.000 


- -C.CC1844 


32.0CC 


-0.001736 


33.CC0 


-0.001637 


34.000 


-0.001547 


35.000 


-0.001463 


36.CC0 


-0.001386 


37.CC0 


-C.C01315 


38.0C0 


-0.001250 


39.000 


-0.001189 


4C.CC0 


-0.001133 


41. COO 


-C. 001080 


42.CCC 


-0.001031 


43.000 


-0.000985 


44.CQC 


-C.0C0943 


45. COO 


-O.C00903 


46. COO 


-0 .OCO 86 5 


47. COO 


-C. 000830 


48.000 


-C. 000797 


49. COO 


-0.000766 


50. COO 


-C.OOC736 
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APPENDIX B 



COMPUTER PROGRAM 



A. ZEROS REGION 

This program will plot the realizability region of general bi- 
quadratic function of the form 

+ 2*SIGMAZ*s + OMEGAZ^ 

Z(s) * 

+ 2*SIGMAP*s + OMEGAP^ 

The input to this program is via punched card. One card is required 
for each plot. The input variables are Sigmap and Omegap specified in 
columns 1 - 20 as two fields of F 10.3 format. 

The output is a plot of allowable region for zeros of Z(s) for a given 
pair of poles. The criterion for the plot is that Z(s) must be positive 
real. 

The output may be either one (or both) of two types of plots depend- 
ing on the needs of the programmers. If a fast solution is desired the 
"CALL DRAW" statements must be removed to enable the program to be run 
on QUICKRUN. The result will be a rough plot on the line printer 
output. 

If TnayTmiim accuracy is desired a CALCOMP plot may be obtained by 
making a regular FORTCLG run with the "CALL DRAW" cards in the deck. 

In this case a line printer plot will also be generated, unless the 
"CALL PLOTP" cards are moved- For the CALCOMP run, "REGION. GO— 10 OK" 
must be specified on the EXEC card. 

The routines DRAW and PLOTP are called by this program. Both are 
in the MPSLIB at NPS. Neither a copy of the decks nor special job 
control language is needed to run this program on the NPS IBM 360. 
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Computer program for the realizability region 
of general biquadratic function. 




i 



The computer program is shown in Figure A- 14. 



B. SENSITIVITY FUNCTION 

This program will plot the sensitivity of general biquadratic 
function with respect to the changing of its components at a particular 
frequency. 

The equation used in this program is. 



(joi.Rp 



W (T(jo),R^ - 0) - T(ju,R^ = INF.) 

(W + R^)^ 



where T(jo),Rp, W, Kju.R^ = 0), T(jo),R 2 = INF.) and R^ correspond to 
Figure 17 in Chapter IV. E. 

The output is a plot of the magnitude of the sensitivity function 
above versus R^* 

The input of this program is via punched card. One card is required 
for each plot. The input variable is frequency in rad/sec. specified 
in columns 1 - 10 as one field of F 10 .*3 format. 

The output is a rough plot on line printer output. 

The routine "CALL PLOT?" is called by this program. 

The computer program is shown in Figure B-11. 
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Computer program for sensitivity function of 
given particular biquadratic function. 
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